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We consider the PT symmetric flat Friedmann model of two scalar fields with positive 
kinetic terms. While the potential of one ("normal") field is taken real, that of the other 
field is complex. We study a complex classical solution of the system of the two Klein- 
Gordon equations together with the Friedmann equation. The solution for the normal 
field is real while the solution for the second field is purely imaginary, realizing classically 
the "phantom" behavior. The energy density and pressure are real and the correspond- 
ing geometry is well-defined. The Lagrangian for the linear perturbations has the correct 
potential signs for both the fields, so that the problem of stability does not arise. The 
background dynamics is determined by an effective action including two real fields one 
normal and one "phantom" . Remarkably, the phantom phase in the cosmological evolu- 
tion is transient and the Big Rip never occurs. Our model is contrasted to well-known 
quintom models, which also include one normal and one phantom fields. 
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1. Introduction 

Complex (non-Hermitian) Hamiltonians with PT symmetry have been vigorously 
investigated in quantum mechanics and quantum field theorjffl. A possibility of ap- 
plications to quantum cosmology has been pointed out in Ref. [21 In the present 
contribution we mainly focus attention on complex field theory. We explore the 
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use of a particular complex scalar field Lagrangian, whose solutions of the classi- 
cal equations of motion provide us with real physical observables and well-defined 
geometric characteristics. 

The interest of our approach is related to its focusing on the intersection between 
two important fields of research: FT symmetric quantum theory and the cosmology 
of dark energy models. 

As is well known the discovery of cosmic acceleratiorP has stimulated an inten- 
sive study of models of dark energjEI responsible for the origin of this phenomenon. 
Dark energy is characterized by a negative pressure whose relation to the energy 
density w = p/e is less than —1/3. Moreover, if this relation happens to be less 
than — 1 such kind of dark energy is called "phantom" dark energji^. The condition 
w > —1 is called the dominant energy condition. The breakdown of this condition 
implies the so called super-acceleration of cosmological evolution, which in some 
models culminates approaching a new type of cosmological singularity called Big 
Rip'^. While the cosmological constant {w = —1) is still a possible candidate for the 
role of dark energy, there are observations giving indications in favor of the models 
where the equation of state parameter w not only changes with time, but would be 
less than — 1 ncP. A standard way to introduce the phantom energy is to consider a 
scalar field with the negative sign of the kinetic energy term, usually called "phan- 
tom" scalar field. However such a model has been reasonably criticized insofar as it 
is unstable with respect to linear perturbations . Constraints on phantom particles 
have been discussed in Ref. [TOl 

In the present paper we propose a cosmological model inspired by FT symmetric 
thcor J- -'^ , choosing potentials so that the equations of motion have classical phantom 
solutions for homogeneous and isotropic universe. Meanwhile quantum fluctuations 
have positive energy density and this ensures the stability around a classical back- 
ground configuration. Thus our study is inspir ed by samples of potentials providing 
the real energy spectrum bounded from belowffEll^IISIH, 

The FT symmetric approach to the extension of quantum physics consists in 
the weakening of the requirement of Hermiticity, while keeping all the physical 
observables real. It was shown that the axioms of quan tum theory are maintained 
if the complex extension preserves (C)FT symmetryi4l. There have also been some 
attempts to apply the FT symmetric formalism to cosmologj02l do not enter 
into the discussion about the equivalence between the FT symmetric theories and 
Hermitian theories with non-local Hamiltonians. The reader is refereed to review 
papers in Refs. HU [151 

We consider the complex extension of matter Lagrangians requiring the reality 
of all the physically measurable quantities and the well-definiteness of geometrical 
characteristics. We would like to underline here, that we work with real space-time 
manifolds. The attempts to use complex manifolds for studying the problem of dark 
energy in cosmology were undertaken in Refs. 16. 

We start with the flat Friedmann model of two scalar fields with positive kinetic 
terms. The potential of the model is additive. One term of the potential is real. 
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while the other is complex and PT symmetric. We find a classical complex solution 
of the system of the two Klein-Gordon equations together with the Friedmann equa- 
tion. The solution for one field is real (the "normal" field) while the solution for the 
other field is purely imaginary, realizing classically the "phantom" behavior. More- 
over, the effective Lagrangian for the linear perturbations has the correct potential 
signs for both the fields, so that the problem of stability does not arise. However, 
the background (homogeneous Friedmann) dynamics is determined by an effective 
action including two real fields one normal and one phantom. As a byproduct, we 
notice that the phantom phase in the cosmological evolution is inevitably transient. 
The number of phantom divide line (PDL) crossings, (i.e. events such that the ratio 
w between pressure and energy density passes through the value —1) can be only 
even and the Big Rip never occurs. The avoidance of Big Rip singularity constitutes 
an essential difference between our mo del a nd well-known quintom models, includ- 
ing one normal and one phantom fieldfl^. The other differences will be discussed 
in more detail later. 

The paper is organized as follows: in Sec. 2 we describe the cosmological model 
we want to analyze, together with a brief explanation of PT symmetric quantum 
mechanics; in Sec. 3 we present the results of the qualitative analysis and of the 
numerical simulations for the dynamical system under consideration; conclusion and 
perspectives are presented in the last section. 

2. Phantom and stability 

We shall study the fiat Friedmann cosmological model described by the metric 

ds^ = dt^ - a^{t)dl^, 

where a(t) is the cosmological radius of the universe. The dynamics of the cos- 
mological evolution is characterized by the Hubble variable 

ft. = -, 
a 

which satisfies the Friedmann equation 

= £, (1) 

where e is the energy density of the matter populating the universe. 

We consider the matter represented by scalar fields with complex potentials. 
Namely, we shall try to find a complex potential possessing the solutions of clas- 
sical equations of motion which guarantee the reality of all observables. Such an 
approach is inspired by the quantum theory of the PT symmetric non-Hermitian 
Hamiltonians, whose spectrum is real and bounded from below. Thus, it is natural 
for us to look for Lagrangians which have consistent counterparts in the quantum 
theory. 

Let's elucidate how the phantom-like classical dynamics arises in such La- 
grangians and for this purpose choose the one-dimensional PT symmetric poten- 
tial of anharmonic oscillator V'''^^^\q) — Xq^{iqy, < e < 2 which has been 
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to possess a real energy spectrum bounded from below. For illus- 
tration, we restrict ourselves with e = 2,V'^^\q) ~ — Ag^. The classical dynamics 
for real coordinates q{t) offers the infinite motion with increasing speed and energy 
or, in the quantum mechanical language, indicates the absence of bound states and 
unboundness of energy from below. However, just there is a more consistent solution 
which, at the quantum level, provides the real discrete energy spectrum, c erta inly, 
bounded from below. It has been proven, first, by means of path Integra and, 
further on, by means of the theory of ordinary differential equation|IS!. In fact, this 
classically "crazy" potential on a curve in the complex coordinate plane generates 
the same energy spectrum as a two-dimensional quantum anharmonic oscillator 
V^^^q) = Mqi + qt)/^ ■^ith real coordinates qi^2 in the sector of zero angular 
momenturrPJJ. Although the superficially unstable anharmonic oscillator is well de- 
fined on the essentially complex coordinate contour any calculation in the style 
of perturbation theory (among them the semiclassical expansion) proceeds along 
the contour with a fixed complex part (corresponding to a "classical" solution) and 
varying unboundly in real direction. In particular, the classical trajectory for V^'^^ {q) 
with keeping real the kinetic, (g)^ and potential, —Xq'^{t) energies (as required by 
its incorporation into a cosmological scenario) can be chosen imaginary, 

q = it e - -2AC^ i^^C- AC^ C > 

which obviously represents a bounded, finite motion with |^| < (C/A)^/^. Such a mo- 
tion supports the quasi-classical treatment of bound states with the help of Bohr's 
quantization. Evidently, the leading, second variation of the Lagrangian around this 
solution, q{t) — i^{t) + Sq{t) gives a positive definite energy, 

£(2) ^ p^t)Sq{t) -H, H= ^pHt) + 12e{t){dqit))^ 

realizing the perturbative stability of this anharmonic oscillator in the vicinity of 
imaginary classical trajectory. It again reflects the existence of positive discrete 
spectrum for this type of anharmonicity. However the classical kinetic energy — 
is negative, i.e. it is phantom-like. 

In a more general. Quantum Field Theory setting let us consider a non-Hermitian 
(complex) Lagrangian of a scalar field 

L^^d^cPdf^r -V{<l>,r), (2) 
with the corresponding action, 

5(^, 0*, 5) = I d'x^/^\{L + ^Rig)), (3) 

where ||g|| stands for the determinant of a metric g'"' and R{g) is the scalar curvature 
term and the Newton gravitational constant is normalized to S/Stt to simplify the 
Friedmann equations further on. 

We employ potentials V{^, $*) satisfying the invariance condition 

(l^($, $*))*= t/($*,$), (4) 
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while the condition 

$*))*= (5) 

is not satisfied. This condition represents a generalized requirement of [C)PT sym- 
metry. 

Let's define two real fields, 

+ x^^i'^-'^*)- (6) 

Then, for example, such a potential can have a form 

$*) = + $*,$-$*) = V{(j), ix), (7) 

where V{x, y) is a real function of its arguments. In the last equation one can 
recognize the link to the so called {C)PT symmetric potentials if to supply the field 
X with a discrete charge or negative parity. When keeping in mind the perturbative 
stability we impose also the requirement for the second variation of the potential 
to be a positive definite matrix which, in general, leads to its PT-symmetrjEl. 

Here, the functions and x appear as the real and the imaginary parts of the 
complex scalar field $, however, in what follows, we shall treat them as independent 
spatially homogeneous variables depending only on the time parameter t and, when 
necessary, admitting the continuation to complex values . 



3. Cosmological solution with classical phantom field 

It appears that among known PT symmetric Hamiltonians (Lagrangians) possessing 
the real spectrum one which is most suitable for our purposes is that with the 
exponential potential. It is connected with the fact that the properties of scalar 
field based cosmological models with exponential potentials are well studied"*^^. In 
particular, the corresponding models have some exact solutions providing a universe 
expanding according to some power law a{t) = oo^'. We shall study the model with 
two scalar fields and the additive potential. Usually, in cosmology the consideration 
of models with two scalar fields (one normal and one phantom, i.e. with the negative 
kinetic term) is motivated by the desire of describing the phenomenon of the so 
called phantom divide line (PDL) crossing. At the moment of the phantom divide 
line crossing the equation of state parameter w = p/e crosses the value w = —1 and 
(equivalently) the Hubble variable h has an extremum. Usually the models using 
two fields are called "quintom models"!^. As a matter of fact the PDL crossing 
phenomenon can be described in the models with one scalar field, provided some 
particular potentials are 

choseiPED 

or in the models with non-minimal coupling 
between the scalar field and gravitjE2l However, the use of two fields make all the 
considerations more simple and natural. In this paper, the necessity of using two 
scalar fields follows from other requirements. We would like to implement a scalar 
field with a complex potential to provide the effective phantom behavior of this field 
on some classical solutions of equations of motion. Simultaneously, we would like to 
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have the standard form of the effective Hamihonian for linear perturbations of this 
field. The combination of these two conditions results in the fact the background 
contribution of both the kinetic and potential term in the energy density, coming 
from this field are negative. To provide the positivity of the total energy density 
which is required by the Friedmann equation ([1]) we need the other normal scalar 
field. Thus, we shall consider the two-field scalar Lagrangian with the complex 
potential 

12 • 2 

L - y + Y - Ae"^ + Be'^x, (8) 

where A and B are real, positive constants. This Lagrangian is the sum of two 
terms. The term representing the scalar field is a standard one, and it can gener- 
ate a power-law cosmological expansiorP^. The kinetic term of the scalar field x is 
also standard, but its potential is complex. In quantum mechanics dynamical sys- 
tems like the latter have been studied in the framework of PT symmetric quantum 
mechanicJi] and, in particular, the exponential potential has been analyzed in great 
detail in Ref. 1131 The most important feature of this potential is that the spectrum 
of the corresponding Hamiltonian is real and bounded from below, provided correct 
boundary conditions are assigned. 

Inspired by this fact we are looking for a classical complex solution of the system, 
including two Klein-Gordon equations for the fields and x- 

^ + 3h^ + Aae"'^ = 0, (9) 
X + 3hx- iBPe'l^^ = 0, (10) 

and the Friedmann equation 

The classical solution which we are looking for should provide the reality and pos- 
itivity of the right-hand side of the Friedmann equation (|lip . The solution where 
the scalar field (j) is real, while the scalar field x is purely imaginary 

X = it e real, (12) 

uniquely satisfies this condition. Moreover, the Lagrangian ([5]) evaluated on this 
solution is real as well. This is remarkable because on homogeneous solutions the 
Lagrangian coincides with the pressure, which indeed should be real. 

Substituting the equation (|12p into the Friedmann equation (jlip we shall have 

h^^^-^l + Ae"^ - Be-^i. (13) 

Hence, effectively we have the Friedmann equation with two fields: one (0) is a 
standard scalar field, the other (^) has the phantom behavior, as we pointed out 
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above. In the next section we shall study the cosmological dynamics of the (effective) 
system, including p3)l . Q and 

i + 3hi~ Bf3e-'^^ ^0. (14) 

The distinguishing feature of such an approach rather than the direct construc- 
tion of phantom Lagrangians becomes clear when one calculates the linear pertur- 
bations around the classical solutions. Indeed the second variation of the action for 
the field x gives the quadratic part of the effective Lagrangian of perturbations: 

Leff = \5x - - Bp^e^^^^'{5x)\ (15) 

where xo is a homogeneous purely imaginary solution of the dynamical system 
under consideration and V stands for spatial gradients . It is easy to see that on 
this solution, the effective Lagrangian (fT5|) will be real and its potential term has 
a sign providing the stability of the background solution with respect to linear 
perturbations as the related Hamiltonian is positive, 

H''eff = \^^^ + \{^{^X)?+Bp^e-K»{5x)\ Sn ^ Sx- (16) 

The Hamiltonian of the metric perturbations should be naturally added to the 
above formulae. This Hamiltonian includes two types of terms: bilinear terms in 
metric perturbations and mixed terms , including both the metric and scalar field 
perturbations. The part of the second-order action, including terms bilinear in the 
metric perturbations can be represented in the following form l -^3 | 24 [ 



C(2) ^ 
9, 9 



+A(r)A'(r) [i(/i^'^ - /i^/io.) + ^(/i''"^,^ - K»\ (17) 

-^{A'\T)+A{T)A"{r))h'>^h], (18) 

where the "prime" means the derivative with respect to conformal time r and the 
background metrics has the form 

while the full metric is 

9tiv{x) = A^{T){ri^^ + h^^{x)). 

The above expression (fT8|) was obtained by the proper integration by parts and 
taking into account the equations of motion B3|24 | r^j^^ only dangerous term here 
is the last "massive" term, which can be removed by the choice of gauge fixing 
condition 



/i"" = 0. 



(19) 
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The second-order action, including mixed terms (the metric perturbations and the 
perturbations of the field x) has the following form: 

= J d'xA'ir) |x' (i/.,o - - (x" + 2^X') h''] <5x, (20) 

where again we have used the equation of motion (the Klein-Gordon equation for the 
field X- This term (pil)) appears annoying because it contains the first derivative of 
the term Be^^^, which is imaginary. However, one can show that by a proper choice 
of the gauge condition this imaginary term can be eliminated. Indeed, the gauge 
condition eliminate from ((20|) the term, proportional to hP^ . Then remains 
the term proportional to x' multiplied by the linear combination of the metric 
perturbations. This combination can be annihilated by the following choice of the 
second gauge condition: 

\h., - h\ = 0. (21) 

This condition is compatible with the first gauge condition p^ . Thus, we have 
shown that the stability with respect to perturbations in our model can be proved 
by fixing a proper gauge condition. 

Let us list the main differences between our model and quintom models, using 
two fields (normal scalar and phantom) and exponential potentialJ^. First, we be- 
gin with two normal (non-phantom) scalar fields, with normal kinetic terms, but 
one of these fields is associated to a complex (PT symmetric) exponential poten- 
tial. Second, the (real) coefficient multiplying this exponential potential is negative. 
Third, the background classical solution of the dynamical system, including two 
Klein-Gordon equations and the Fricdmann equation, is such that the second field 
is purely imaginary, while all the geometric characteristics are well-defined. Fourth, 
the interplay between transition to the purely imaginary solution of the equation 
for the field x s-^id the negative sign of the corresponding potential provides us with 
the effective Lagrangian for the linear perturbations of this field which have correct 
sign for both the kinetic and potential terms: in such a way the problem of stability 
of the our effective phantom field is resolved. Fifth, the qualitative analysis of the 
corresponding differential equations, shows that in contrast to the quintom models 
in our model the Big Rip never occurs. The numerical calculations confirm this 
statement. 

In the next section we shall describe the cosmological solutions for our system 
of equations. 

4. Cosmological evolution 

First of all notice that our dynamical system permits the existence of cosmologi- 
cal trajectories which cross PDL. Indeed, the crossing point is such that the time 
derivative of the Hubble parameter 

h^-\{^''-e) (22) 
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is equal to zero. We always can choose = &i t = tpoL provided the values of 
the fields (j){tpDL) and (,{tpDL) are chosen in such a way, that the general potential 
energy Ae°"^ — Be''^^ is non- negative. Obviously, tpoL is the moment of PDL 
crossing. However, the event of the PDL crossing cannot happen only once. Indeed, 
the fact that the universe has crossed phantom divide line means that it was in 
effectively phantom state before or after such an event, i.e. the effective phantom 
field ^ dominated over the normal field (j). However, if this dominance lasts for a long 
time it implies that non only the kinetic term — ^^/2 dominates over the kinetic term 
(/)^/2 but also the potential term — _Bexp(— should dominate over Aexp(a0); 
but it is impossible, because contradicts to the Friedmann equation (fT3|) . Hence, 
the period of the phantom dominance should finish and one shall have another 
point of PDL crossing. Generally speaking, only the regimes with even number of 
PDL crossing events are possible. Numerically, we have found only the cosmological 
trajectories with the double PDL crossing. Naturally, the trajectories which do 
not experience PDL crossing at all also exist and correspond to the permanent 
domination of the normal scalar field. Thus, in this picture, there is no place for 
the Big Rip singularity as well, because such a singularity is connected with the 
drastically dominant behavior of the effective phantom field, which is impossible 
as was explained above. The impossibility of approaching the Big Rip singularity 
can be argued in a more rigorous way as follows. Approaching the Big Rip, one has 
a growing behavior of the scale factor a{t) of the type a{t) ^ {tBR — i) where 
q > 0. Then the Hubble parameter is 

h{t) = (23) 

and its time derivative 

h{t) 



{tBR-ty 

Then, according to Eq. (HH), 

Substituting Eqs. (|24p . into the Friedmann equation p^ . we come to 

q q 



{tBR - tf i{tBR - tf 



In order for this to be satisfied and consistent, the potential of the scalar field 
(j) should behave as l/{tBR — t)"^- Hence the field should be 

(^ = (j)a--\n{tBR-t), (25) 
a 

where 0o is an arbitrary constant. Now substituting Eqs. and ([^5]) into the 
Klein-Gordon equation for the scalar field ((> (l9|), the condition of the cancelation 
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of the most singular terms in this equation which are proportional to \/{tBR — t)'^ 
reads 



This condition cannot be satisfied because all the terms in the left-hand side of 
Eq. (j26p are positive. This contradiction demonstrates that it is impossible to reach 
the Big Rip. 

Now we would like to describe briefly some examples of cosmologies contained 
in our model, deduced by numerical analysis of the system of equations of motion. 
In Fig. [T]a double crossing of PDL is present. The evolution starts from a Big Bang- 
type singularity and goes through a transient phase of super-accelerated expansion 
("phantom era"), which lies between two crossings of PDL. Then the universe un- 
dergoes an endless expansion. In the right plot we present the time evolution of the 
total energy density and of its partial contributions due to the two fields, given by 
the equations 



e = Ecj, + s^, 

which clarify the roles of the two fields in driving the cosmological evolution. 
The evolution presented in Fig. [2] starts with a contraction in the infinitely remote 
past. Then the contraction becomes superdecelerated and turns later in a super- 
accelerated expansion . With the second PDL crossing the "phantom era" ends; 
the decelerated expansion continues till the universe begins contracting. After a 
finite time a Big Crunch-type singularity is encountered. From the right plot we 
can clearly see that the "phantom era" is indeed characterized by a bump in the 
(negative) energy density of the phantom field. 

In Fig.[3]the cosmological evolution again begins with a contraction in the infinitely 
remote past. Then the universe crosses PDL: the contraction becomes superdecel- 
erated until the universe stops and starts expanding. Then the "phantom era" ends 
and the expansion is endless. 

In Fig. |4] the evolution from a Big Bang- type singularity to an eternal expansion is 
shown. The phantom phase is absent. Indeed the phantom energy density is almost 
zero everywhere. 

4.1. The cosmological dynamics in the presence of dust like matter 

Trying to make our model more realistic we can try to describe which changes 
does it undergo in the presence of some matter. The natural candidate for the role 
of matter is a dust-like perfect fluid with the vanishing pressure. In this case the 
Friedmann equation (|13p acquires the form 



2 + 6q + Ac? exp(a(/)o) = 0. 



(26) 




2 



(27) 
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Fig. 1. (left) Plot of the Hubble parameter representing the cosmological evolution. The evolution 
starts from a Big Bang-type singularity and goes through a transient phase of superaccelerated 
expansion ("phantom era"), which lies between two crossings of PDL (when the derivative of h 
crosses zero). Then the universe expands infinitely, (right) Plots of the total energy density (blue), 
and of the energy density of the normal field (purple) and of the phantom one (green). 




Fig. 2. (left) The evolution starts with a contraction in the infinitely remote past. At the point first 
PDL crossing the contraction becomes superdecelerated and turns in a superaccelerated expansion 
when h crosses zero. The second PDL crossing ends the "phantom era" ; the decelerated expansion 
continues till the universe begins contracting. In a finite time a Big Crunch-type singularity is 
reached, (right) Plots of the total energy density (blue), and of the energy density of the normal 
field (purple) and of the phantom one (green) . 

while the second Friedmann equation (|22p becomes 

where the constant Co characterizes the quantity of the dust-Uke matter in a universe 
under consideration. 

Obviously at the early stage of the cosmological evolution, this dust contribution 
into the Friedmann equations dominates over dark energy contribution coming from 
our scalar fields. So, at early stage of the evolution (close to the Big Bang epoch) the 
effects like phantomization of the matter content of the universe and phantom divide 
line crossing are impossible. They are postponed until the time when the expansion 
of the universe is so large that this additional term becomes insignificant. However, 
one can say that its presence even underlines some peculiar effects of the considered 
model. Indeed, it makes stronger the non-phantom part of the matter content of 
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Fig. 3. (left) The cosmological evolution begins with a contraction in the infinitely remote past, 
with the first PDL crossing the contraction becomes superdecelerated until the universe stops 
{h = 0) and starts expanding. With the second crossing the "phantom era" ends and the expansion 
continues infinitely, (right) Plots of the total energy density (blue), and of the energy density of 
the normal field (purple) and of the phantom one (green). 




J — , — , — , — , — , — , — , — \ — , — , — , — , — , — , — , — , — , — ^ |. J — , — , — , — , — - - I , t 

-0.004 -0.002 0.002 0.004 -0.004 -0.002 0.002 0.004 

Fig. 4. (left) Evolution from a Big Bang-type singularity to an infinite expansion, without any 
crossing of PDL. This evolution is thus guided by the "normal" field <p. (right) Plots of the total 
energy density (blue), and of the energy density of the normal field (purple) and of the phantom 
one (green). Notice that the energy density of the phantom field (green) is very close to zero, thus 
the total energy density is mainly due to the standard field. 

the universe and hence return to the normal non-phantom case occurs more quickly. 
Naturally, the Big Rip singularity is avoided in this model. 

5. Conclusion 

As was already said above the data are compatible with the presence of the phantom 
energy, which, in turn, can be in a most natural way realized by the phantom scalar 
field with a negative kinetic term. However, such a field suffers from the instability 
problem, which ma kes it vulnerable. Inspired by the development of FT symmetric 
quantum theorjEI -^yg introduced the FT symmetric two-field cosmological model 
where both the kinetic terms are positive, but the potential of one of the fields is 
complex. We studied a classical background solution of two Klein-Gordon equations 
together with the Friedmann equation, when one of this fields (normal) is real while 
the other is purely imaginary. The scale factor in this case is real and positive just 
like the energy density and the pressure. The background dynamics of the universe 
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is determined by two effective fields - one normal and one phantom, while the 
Lagraiigian of the linear perturbations has the correct sign of the mass term. Thus, 
so to speak, the quantum normal theory is compatible with the classical phantom 
dynamics and the problem of instability is absent. 

As a byproduct of the structure of the model, at variance with quintom models, 
the phantom dominance era is transient, the number of the phantom divide line 
crossings is even and the Big Rip singularity is excluded. In all our numerical sim- 
ulations we have seen only two phantom divide line crossings. However, one cannot 
exclude that at some special initial conditions, or, in the models of the kind con- 
sidered in the present paper, but with more complicated potentials, it is possible to 
see more than two phantom divide line crossings. It would be especially interesting 
to find the potentials or/and initial conditions providing an oscillatorial behaviour 
with infinite phantom divide line crossings. 
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